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• On the basis of the kinetic equation with selfconsistent relativistic mean fields acting on baryons 

, and antibaryons, we study dynamical response of the nucleus to an antiproton implanted in its 



interior. By solving numerically the time-dependent Vlasov equation, we show that the compressed 
state is formed on a rather short time scale of about 4-^ 10 fm/c. This justifies the assumption, that 
the antiproton annihilation may happen in the compressed nuclear environment. The evolution 
of the nucleus after antiproton annihilation is described by the same kinetic equation including 

(N 

' collision terms. We show, that nucleon kinetic energy spectra and the total invariant mass distri- 

^ : 

' butions of produced mesons are quite sensitive observables to the antiproton annihilation in the 

compressed nucleus. 
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I. INTRODUCTION 

As has been shown recently in Refs. , an antiproton implanted in a heavy nucleus 
serves as an attractor for surrounding nucleons that can lead to a sizable increase of the 
central nucleon density. This effect is caused by the strong attractive scalar and vector 
potentials acting on the antiproton, as follows from the G-parity transformation of nuclear 
potentials Correspondingly, the antiproton also creates an attractive potential acting on 
nucleons. This leads to the concentration of nucleons around the antiproton and, as result, 
to a considerable increase of the nucleon density. 

Within the relativistic mean field (RMF) model, the G-parity transformed nuclear optical 
potential is about —700 MeV at the normal nuclear matter density po = 0.148 fm~^, while 
a phenomenological value of an antiproton optical potential is limited within the range of 
-(100 - 350) MeV Q, H, fl Therefore, in order to fit the empirical optical potential, 
the antiproton coupling constants with a-, uj- and p-meson fields should be reduced with 
respect to the values given by the G-parity transformation. The RMF calculations with 
reduced coupling constants 2] still show quite strong compressional effects for light and 
medium nuclei. 

An important question, which arises here is whether the compression process is fast 
enough to develop before the p-annihilation. The total pp-annihilation cross section in 
vacuum can be parameterized at low relative velocities frei as 

<n=C^+— , (1) 
frel 

n 

where G = 38 mb and D = ?>h mb-c [9,]. Using these numbers we can estimate the life time 
of an antiproton inside the nuclear matter at normal density: 

Tann ^ ^ ^ 2 fm/c . (2) 

Po C" ann f rel 

This is, of course, a very short time in nuclear scale. However, as argued in Ref. j^, this 
time can become much longer, up to 20 fm/c, for deeply bound antiprotons due to the phase 
space suppression factors. Therefore, the compression effects can, in-principle, show up in 
p-nuclear interactions. 

In the present work, we apply a dynamical transport model in order to study the formation 
and decay of the compressed p-nuclear system. Our calculations are based on the Giessen 
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Boltzmann-Uehling-Uhlenbeck (GiBUU) model lOj, which has been recently supplemented 



by the relativistic mean fields 



ll| . Apart from collision terms, the GiBUU model solves the 



coupled (through the mean fields) Vlasov equations for nucleon and antiproton phase space 



distribution functions. As well known 



12| , the Vlasov equation provides a semiclassical limit 



of the time-dependent Hartree-Fock calculations. Thus, the compressional effects found in 
Refs. [1,2] should also be reproduced as a static solution of the coupled Vlasov equations. 

It will be demonstrated that the compression process is characterized by the time scale 
which is comparable with the p life time in nuclear environment. Thus, p has, indeed, a 
chance to annihilate inside the compressed nucleus. We will show, that the p-annihilation 
in a compressed nucleus should lead to the collective expansion of the residual nuclear 
system. The appearance of the high-energy tails in the kinetic energy spectra of the emitted 
nucleons is predicted. The distributions in the total invariant mass of produced mesons 
reveal a noticeable shift toward lower invariant masses, when the annihilation takes place 
inside the compressed nucleus. 

The annihilation of slow antiprotons inside heavy nuclei was, first, proposed by Rafelski 



13| as a unique opportunity to study nuclear matter in unusual conditions. Later, Cahay 



et al. Im studied the p annihilation inside nuclei within an intranuclear cascade model. In 
Ref. [^1, antiproton annihilation events into pions at the center of '^°Ca and ^°®Ag nuclei 
were simulated. The mean field effects were, however, completely neglected in 14 1. 

In Sect. II, we describe the theoretical model applied in calculations. Sect. Ill contains 
the results of the time evolution study for the compression and explosion dynamics. In Sect. 
IV, we propose several observable signals sensitive to the p-annihilation in the compressed 
nucleus. The summary and outlook are given in Sect. V. 



II. THE MODEL 



In calculations, we apply the GiBUU model developed in Giessen Universit y. For the 
detailed description and related references, we refer the reader to the web page [lO|, where 
the new version of the model is presented. Below, we mostly describe the new features 
implemented in the present work. 
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A. Relativistic mean fields 



Below we consider a system composed of an antinucleon interacting with baryons. This 
system is described by the RMF Lagrangian of the following form [2], [l^ : 

j=B,N 

+ ^d^ad^a - U{a) - \f,,F^^ + \mlu,u:^ , (3) 

where ipj are the baryon {j = B ^ N, N*, A, Y) and antinucleon (j = N) fields, respec- 
tively; a is the isoscalar-scalar meson field (J*^ = 0"'", J'^ = 0"'"); cu^ is the isoscalar- vector 
meson field {I^ = 0^, J'^ = and F^,^ = S^cj,^ — d^uj^. Here A^* and A denotes, re- 
spectively, the isospin 1/2 and 3/2 nonstrange baryonic resonances, and Y stands for the 



S = —1 baryons explicitly propagated in the GiBUU model [10|. In the case of the spin 3/2, 
5/2 and 7/2 baryonic resonances, their fields ipj carry also one or more vector indices, which 
are dropped in Eq.(l3]) and below for brevity. When appropriate, the covariant summation 
is assumed over these indices. For simplicity, the isovector and electromagnetic terms are 
disregarded in ([3]). The selfinteractions of the cr-field are included in ([3]) via the term Uja) 



in order to avoid an unrealistically high compressibility coefficient of the nuclear matter 



16|: 



Uia) = \mla' + ^g,a' + \g,a^ . (4) 

Some comments are in order to gain more insight into Eq. Following Ref. j^, 

the antinucleon field ipfi in the Lagrangian density ([3]) is represented in terms of wave 
functions of physical antinucleons. These wave functions can be obtained by the G-parity 
transformation acting on the wave functions of the Dirac sea nucleons (see Ref. [3| for 
details), which appear in the relativistic description of the nucleon 17 1. By applying the 
same transformation, the nonlinear RMF Lagrangian of Refs. 15|, |l6| (neglecting terms 
responsible for the baryon-antibaryon annihilation) can be expressed as (E]) with the following 
relations between coupling constants: 

9u)N = —QuiNi QaN = 9aN ■ (5) 

The relations ([S]) are satisfied if the physical system would be exactly symmetric with respect 
to the G-parity transformation. However, this is not necessary to be true in a many-body 



system (2 



, |3| . The reason is that the concept of the G-parity symmetry is strictly applicable 



on the level of the elementary processes only. However, the RMF Lagrangian ([3]) is dealing 
with the effective interactions, which are usually tuned to describe the bulk properties of 
the nuclear medium and/or the properties of some selected nuclei. Due to the many-body 
effects, such as the Pauli blocking or mixed scalar-vector terms in the scattering amplitudes, 
these effective interactions may not obey the exact G-parity symmetry anymore. To take 
into account possible deviations from the G-parity symmetry, we introduce an overall scaling 
of the antinucleon-meson coupling constants with respect to the values given by ([5]) (see Ref. 

duiN = -iduN-, 9aN = idaN , (6) 

where < ^ < 1 is a scaling factor. 

Throughout the paper, we consider two options for the scaling factor of the antinucleon- 
meson coupling constants: = 1, motivated by the G-parity, and C, = 0.3, which is in 
a better agreement with the empirical pA optical potential. For other baryonic fields we 
put in the present work, for simplicity, the same coupling constants as for the nucleon : 

9ujN* = 9ujA = 9ujY = 9ujN, 9crN* = 9(7 A = 9crY = 9(tN- , , 

All calculations have been performed emloying the NL3 parameterization [15] of the 
RMF model. This parameterization provides quite reasonable nuclear matter properties: 
the binding energy 16.299 MeV/nucleon, the compressibility coefficient K = 271.76 MeV 
and the nucleon effective mass = 0.60 m^r at po- Moreover, the NL3 parameterization 
reproduces the ground state properties of spherical and deformed nuclei very well [15|. 

The Dirac equations of motion for baryons have the following form: 

{-f'^iid^ - 9^jU^) - m))iljj = , (7) 

where 

m* = nij + g^ja (8) 

is the effective (Dirac) mass. 

Within the mean field approximation the a- and tu-fields are treated classically. They 
satisfy the (nonlinear) Klein-Gordon-like equations with the source terms due to coupling 
to baryons and an antinucleon: 

dud^^ + ^ = - E 9.,Ps^ , (9) 

j=B,N 

[d.d^ + ml) u;^= Y: a.jJbj ' (10) 



where psj =< ^Jj^Jj > is the partial scalar density and jj^- =< ipjl'^i'j > is the partial baryon 
current. Equation (|T0|) has to be supplemented by the four-transversality condition 

d^oo'' = . (11) 



B. Covariant kinetic equations 



Instead of solving the Dirac equations ([7]), we will describe the 
dynamics by the coupled set of the semiclassical kinetic equations 



mryons and antinucleon 
1^ 
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2H: 



1 

Po 



P 



-kfl 



d_ 



+ (9^jPlF''^ + 



dm* 



d 



f,{x,pn = iAfB,L 



Ml 



(12) 



dxk J dp*'' 

where = 1,2,3; /x = 0,1,2,3; x = (t,r); and fj{x,p*) is the distribution function (DF) 
in a six- dimensional phase space (r, p*) with being the spatial components of the kinetic 
four- momentum 



p*>' = p'' - g^,uj>' 



(13) 



The baryons and antinucleon are assumed to be on the respective effective mass shells: 



P 



*\2 



(14) 



The l.h.s. of Eq. (fT2l) describes the propagation of the j-th type particles in the classical a- 
and tu-fields. The r.h.s. of Eq. f ll2l) is a collision integral, which represents the (in)elastic two- 
body collisions with corresponding vacuum cross sections as well as the resonance decays. 
The complete description of the collision integral structure, in-particular, the differential 
elementary cross sections included into the GiBUU model can be found in 10|, lll| and in 
refs. therein. The in-medium modification of the baryon-baryon and baryon-meson cross 
sections is neglected in the present work. 

We will apply the full kinetic equations, including collision terms, only to describe the 
post-annihilation evolution of a system. By this reason, the antiproton DF is excluded 
from the collision integral. Instead, we enforce p to annihilate into mesons at some pre- 
selected time (see Sect. HE). Thus, in the present work the collision integral includes 
the nucleon, A(1232) and higher baryon resonances up to the mass of 2 GeV, which can 
be excited in the meson-baryon and baryon-baryon collisions. A possible hyperon forma- 
tion in the processes vrA^ YK and KN —>■ nY is included too. The "valence mesons" 
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M = n, 7], p, a, u, r]', 0, r]c, J/'ip, K, K, K* , .ft'* are explicitly taken into account. They 
are assumed to propagate freely between collisions, i.e. we neglect the mean field potentials 
acting on these mesons. 

The scalar density and the baryon current of the j-th type baryons are expressed in terms 
of DF as follows: 



where gj is the spin-isospin degeneracy factor {g^ = (7^ = 4 , (7a = 16 etc. 



One can show 



21 



221 ] that the kinetic equations (IT^ with the a- and u;-fields evolving 



according to Eqs. fP). ffT(l lead to the continuity equations 

E ^fb, = , d,fm = (17) 
and the energy-momentum conservation 

duT^''' = , (18) 
where the energy-momentum tensor is written as 

,=i:^,M(2^) P 

-g^- i^d^ad^a - U{a) - hxco^d^co^ + ^m>2^ . (19) 

Here we have also included possible contributions of the "valence" mesons M, which can be 
produced at the annihilation. It is assumed that p* = p for the valence mesons. 



Numerical realization 



In order to solve Eq. (1121) numerically, DF is represented by the set of point-like test 
particles: 

/,(x, P^) = ^ E ^(r - r.(t))5(p^ - pK^)) , (20) 

where Nj is the number of physical particles of the type j and n is the number of test 
particles per physical particle (the same for all types j). The test particle positions r, and 



kinetic momenta p* are evolving in time according to the following equations: 

with k = 1,2,3 and yU = 0, 1, 2, 3. It is easy to check that DF fl20|) with and satisfying 
Eqs. ( I2T1) ■ (|22l) gives a formal solution of the Vlasov equation in the case when the collision 
integral in ( |T2l) is equal to zero. Equations ( 12T|) . ( l22i) are equivalent to the Hamiltonian 
equations of motion for the test particle positions and canonical momenta pj: 



(23) 
(24) 



where p° = g^jUJ^ + ^(Pi)^ + (m^)^ is the single-particle energy (see 2l|, l23|). However, it 
is more convenient to propagate in time the test particle kinetic momenta rather than the 
canonical ones, since then Eq. ([9]) for the cr-field decouples from Eq. f fTOj) for the cu-field. 

When the collision integral in Eq. (fT2!) is taken into account, the test particles are 
propagated between the two-body collisions using Eqs. fl?I]) . fl22]) . All calculations have been 
performed in the parallel ensemble mode. In this mode, the two-body collisions are permitted 
between the test particles belonging to the same parallel ensemble only, while the mean field 
is averaged over n parallel ensembles of the test particles propagated simultaneously (see 
Eq. (!20!) ). Therefore, a single parallel ensemble can be considered as a physical event. 

In actual calculations, we have neglected the time derivatives of the meson fields in Eqs. 
(P). (ITU]) . However, the spatial derivatives were treated without any simplifying assumptions. 
The reason for such a strategy is that we are dealing with nuclear systems which have large 
density gradients, but evolving slowly, as compared with the spatial and temporal scales 
involved in the mesonic equations of motion. Indeed, including the temporal gradients would 
lead to the frequent oscillations of the mesonic fields with a period of less than l-n jm, where 
m is the meson mass. This gives the period of 2.5 fm/c (1.5 fm/c) for the a- (cu-) field. 
By taking into account the finite wave lenghts of these oscillations would further reduce 
the periods. The treatment of such oscillations would strongly complicate the numerical 
calculations, in-particular, due to the classical meson field radiation. On the other hand, 
the characteristic periods of the oscillations are significantly smaller than the characteristic 

8 



compression times (4 10 fm/c, see Sect. Ill A below). Therefore, one can approximately 
average-out the mesonic fields with respect to these oscillations, that is actually assumed in 
our model. The a- and a;-fields are, therefore, calculated from the equations 

j=B,N 

{-A + ml)uj^= Y: 9u.,jI:,. (26) 

j=B,N 

Within the same approximation, the energy-momentum tensor has the following form: 
= E TTTTs / ^pVV,(x, P^) + (9^a9V - d>^u;'d^co,)il - 5.o) 

j=B,N,M '^^^^ P 

-9"" (-^(V^)' - U{a) + iv^.Va;^ + ^m^^) . (27) 

The factor (1 — 6^0) in Eq. (1271) refiects the fact, that due to the omission of the time 
derivatives of the meson fields in the Lagrangian density, only the first term in the r.h.s. 
contributes to the three-momentum density {a = 1, 2, 3). 

Although Eqs. fl251) . fl2Bl) are not covariant, they provide a better description of the nuclear 



surface than pure local fields . This improves the stability of a nuclear ground state and 
is more appropriate for studying nuclear response to external hadronic and electromagnetic 
probes. 

Equations (p5l) . (!26|) have been solved numerically by a pply ing the alternating direction 
implicit iterative method of Douglas described in Ref. 2J]. Due to the scalar density 
dependence on the effective mass (see Eqs. (15]). (ITS]) ), additional iterations are needed to solve 
Eq. (l2S]) . In other words, the scalar density has to be computed selfconsistently. To evaluate 
the meson fields, we used a uniform grid in coordinate space with steps A^; = Ay = A^. 
For the systems p^^O, p^'^Ca and p^^^Ph considered below, the grid covered a cubic volume 
with the side of 10, 20 and 30 fm, respectively, centered at the center-of-mass (cm.) of a pA 
system. By numerical reasons, the (5-functions in coordinate space, introduced in Eq. (12Uj) . 
have been replaced by the Gaussians of the width L: 

The width of the Gaussian and the grid step sizes are pure numerical parameters which 
should resolve the coordinate space nonuniformities of the system. In our case the charac- 
teristic space scale is given by the radius of the smallest considered nucleus ^^O, i.e. ^ 3 



fm. On the other hand, in order to have smooth density distributions, the number n of 
test particles per physical particle (see Eq. fl20l) ) should be correlated to the width of the 
Gaussian as n oc L^^. This puts a restriction on too small width due to CPU time increase. 
As an optimum choice, we fixed in the present work A^. = A.^ = = L = 0.5 fm. We 
have realised, however, that there is a rather moderate tendency of increasing maximum 
compression (see discussion in Sect. Ill A below) with decreasing Gaussian width. The 
number of test particles per nucleon was set to n = 1500 in the most of calculations. 

The equations of motion fl2T!) . (!22|) have been solved by applying the second-order in time 
predictor-corrector method ll| with the time step of 0.1 fm/c. This value is small enough 
to resolve the time scale of a few fm/c for the compression processes (see Figs. E] and E] 
below). We have checked, that taking smaller time step does not influence the results. The 
full numerical scheme conserves the total energy with the accuracy of about 5% of the initial 
total binding energy of the pA system. 

D. Initialization 

The nucleons were distributed in coordinate space according to the Woods-Saxon density 
profile. The momenta of nucleons were sampled according to the local Fermi distribution. 
The initial antiproton DF was chosen as a Gaussian wave packet in coordinate and mo- 



mentum space [25|, l26| located at the center of a nucleus {x = y = z = 0): 



Mt = 0, r, p*) = ^ exp{-rV(2a,2) - 2a',p^'} , (29) 



where is the width in coordinate space. Equation (!29i) implies, that the antiproton is at 
rest. The width of the initial antiproton distribution in momentum space is (2cTj.)~^, which 
follows from the uncertainty relation. If not mentioned explicitly, the calculation is done 
with the choice cr^ = 1 fm. This value agrees with results of the static RMF calculations of 
Ref. As for nucleons, the antiproton DF ( !29l) is projected onto test particles according 
to Eq. ( l20l) with the 5-functions in coordinate space replaced by Gaussians. To avoid mis- 
understanding, we note that one should distinguish the width cr,. of the physical antiproton 
spatial distribution in Eq. fl29l) and the width of the test particle Gaussian. 
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E. Propagation and annihilation 



After the initialization, the system of nucleons and antiproton was propagated in time 
according to Eqs. fl^ . fl221) . The meson fields have been calculated by Eqs. fl2S]) . fl2^ with 
the source terms given by the scalar densities fllSp and the baryon currents (ITB]) . In such 
a way, the evolution of the system toward compressed state has been followed. In this 
calculation, the collision term in the r.h.s. of the kinetic equation f|T2|) has been set to zero, 
i.e. we considered a pure mean-field Vlasov dynamics. This was done to see most clearly the 
role of the mean fields. An introduction of the NN and A^A^ elastic collisions would mainly 
lead to a dissipation of the collective energy into heat. As pointed out in Ref. 21], this effect 
is rather small and, therefore, can not change significantly the compression dynamics. 

The reason is, that the elastic collisions are not frequent on the time scale of compression 
(see Figs. [3], H] and [5]) . Indeed, the mean time Tcou between nucleon-nucleon collisions can be 



estimated as Tcoh 



section (c.f. Refs. [271, 



1/ 



[rnO'nn'Vf), where a^N ~ 40 mb is the elastic nucleon-nucleon cross 



28| ) and vp — 0.3c is the Fermi velocity. This gives Tcoh = 3^6 fm/c 



for the nucleon density piv = 2 4- Ipo- The Pauli blocking effect will further increase Tcou. A 
similar estimate can also be done for A^A^ elastic collisions. 

At certain time moment tann, which is an external parameter to our model, we simulated 
the annihilation of an antiproton. This implies that annihilation occurs instantaneously, as 
a single quantum mechanical transition, in distinction to description of this process via the 
collision term in a kinetic equation. In the last case, the antiproton distribution function 
would gradually disappear on the way to the compressed state. The purpose of the present 
work is to look at the strongest possible effect of the nuclear compression on observables. 
Therefore, we let the compressed system to be formed, and simulate the sudden annihilation 
afterwards. The ambiguity in the in-medium annihilation cross sections is taken into account 
by varying the parameter tann- 
in the actual calculations, the annihilation was simulated as follows: For each antipro- 
ton test particle, the closest in coordinate space nucleon test particle was chosen to be the 
annihilation partner. At large enough values of the total in-medium cm. energy ^/s of the 
annihilating pN pair (see below), the annihilation event of the test particle pair into mesons 
was simulated using ,ua..U .node, Q, wMcK Ms been already „np>e,nen.ed ,n the 
GiBUU model lOj earlier. A quark and an antiquark with the same flavour are assumed 
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to annihilate and transfer their total four- momentum to the remaining (anti) quarks. The 
remaining four (anti)quarks form two orthogonal qq jets with equal energies in the cm. 
frame. The jets were hadronized via the Lund string fragmentation model 3l|] in the JET- 
SET version included into the PYTHIA 6.225 program package. The applied annihilation 
model corresponds to the R2 type diagram in classification of Ref. [^], i.e. to the quark 
rearrangement with one qq annihilation vertex. In this sense, the model has some similarity 



3d, 



33j. To illustrate how the model works 



with the two-meson doorway models of Refs. 
we have performed simulations of the pp annihilation. 

Fig. [1] shows the pion multiplicity distribution for the pp annihilation at rest in vacuum 
compared to the data compilation from Refs. [9|, 13J]. The calculated distribution is some- 
what shifted to smaller pion multiplicities with respect to the data: The calculated average 
pion multiplicity < ra^ >^ 4.5 compared with the experimental value of ~ 5.0. We would 
like to remark, that non-vanishing contribution of the = 2 channel in calculations is 
completely due to the final states with other particles: n-nrj (78%), rniKK (14%), nnriri 
(6%) and vrvr + photons (2% before r] decay). In calculations, we took into account t] decays 
into 27 or into final states with pions and disregarded photons afterwards. However, it is 
not clear to us how photons were counted in the data (see also Ref. [32]). 

Fig. [2] shows the calculated charged pion momentum distributions in the cm. frame of 
the annihilating pp pair at rest in vacuum. From the partial contributions of the channels 
with various pion multiplicities we observe, as expected, that the hard(soft) part of the total 
momentum distribution is populated mainly by the low(high) pion multiplicity events. The 
experimental data are described reasonably well, except for the momenta 0.5 < k < 0.7 
GeV/c, where the calculations significantly overestimate the data. 

We believe that the accuracy of the model in describing the data in Figs. [1] and [2] is 
sufficient for the exploratory studies of global observables in the present work. Certainly, 
the improvement of the annihilation model is needed to perform more detailed study of the 
mesonic final states in the annihilation. Below we concentrate more on the in-medium effects 
on the annihilation. 

Due to the mean field, the invariant energy of the annihilating pN pair can be substan- 
tially below the vacuum threshold value of 2m7v. This makes the direct application of the 
JETSET model for the p-annihilation in nuclei physically and numerically problematic. To 
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overcome this difficulty, we introduced the corrected invariant energy as follows 1 ll. |35[| : 



V^corr = - 2(m^ - rriN) , (30) 

where s* = (p| + p^)^. The quantity ^/s^ is a vacuum analog of the total in-medium 
invariant cm. energy with s = {pp + PnY- Provided that ^/s > 4m7r, we have used 
^Scorr ill the JETSET simulation in order to produce the mesonic final states. This lower 
limit of ^/s is due to the fact that the JETSET model does not generate enough direct 27r 
and 37r annihilation final states. 

In order to take into account the in-medium effects, in-particular, to ensure the correct 
in-medium threshold condition y/s > mi + m2 + ... + rrin^^^, where mi, m2, ...,m„^^^ are 
the vacuum masses of the produced mesons, the annihilation event was accepted with the 
probability 

where 



X <5(^)(pp + pN-h-h-...- K^J (32) 

is the invariant phase space volume, ki = {uJi, ki) are the four-momenta of the produced 
mesons satisfying the vacuum mass shell conditions mf = kf, i = 1,2, rimes- Finally, the 
three-momenta of the produced mesons in the cm. frame of the annihilating pN pair were 
multiplied by the common factor adjusted to get the correct in-medium total cm. energy 

In-fact, the way we simulate the in-medium effects Eg. (1311) implies using the vacuum 
matrix elements of the annihilation channels, which are given by the JETSET model, while 
taking into account the in-medium effects in the phase space factors only. Similar procedures 



have been applied earlier in Refs. I2I, [ill, l35 |. 



At 2m7r < < 4m7r, the final 27r or Sir channel was chosen by Monte-Carlo according 
to the probability ratio 

^27r _ ^ $2(v^; m^, m^)$3(2mAr; m^, m^, m^) 
V^n °$2(2miv;m^,m^)$3(v/i;m^,m^,m^) 

where Rr, = 0.152 is the ratio of the 27r and final state probabilities for the pp annihilation 

n 

at rest (see Table VI in Ref. ||). For the zero total charge Q of the annihilating pN pair, the 
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charge states of the outgoing pions were also determined from the data compilation of Ref. 
[2I. Since for Q = ±1 the data are absent, the charges of the Sir final states were determined 
by assuming that the tt^tt^tt^ and tt'^tt+tt^ final channels have equal probabilities. The 
momenta of the outgoing pions were distributed microcanonically according to the available 
two- or three-body phase space. 

After the annihilation is simulated, the residual nucleons and produced mesons were 
propagated in time according to the full kinetic equations f|T2l) . including both the baryonic 
mean fields and collision integrals. This takes into account the entropy production caused 
by the two-body collisions at the expansion stage. Moreover, important processes of the 
meson rescattering and absorption, e.g. vrA^ ^ A — nN or ttN A, AA^ — A^A^ are 
included in the collision integral. These processes influence the observed particle spectra. 

III. TIME EVOLUTION OF BOUND p-NUCLEAR SYSTEMS 
A. Initial compression stage 



As demonstrated in Refs. [1, 2!| by static RMF calculations, a deeply-bound antiproton- 
nucleus system can be significantly compressed as compared with a normal nucleus. Now 
we want to study the real dynamics of such a system starting from the unperturbed nuclear 
ground state at t = 0. 

Fig. [3] (top panels) shows the nucleon and antiproton density profiles calculated at dif- 
ferent times along the axis z drawn through the center of the p^°Ca system. Fig. [3] (bottom 
panels) also shows the nucleon and antiproton potentials Uj = g^ijOJ^ + QajO', j = N, p, along 
the same axis. Left and right panels present results for ^ = 0.3 and ^ = 1, respectively. We 
see that the initial configuration is unstable and the system starts to shrink. Both nucleon 
and antiproton central densities grow quite fast, reaching their maxima within several fm/c. 
In the course of the compression process, the nucleon potential becomes deeper in the case 
of the reduced antiproton coupling constants {C, = 0.3) and does not, practically, change in 
the G-parity motivated case (^ = 1). The antiproton potential deepens quite strongly with 
time for the both sets of the antiproton coupling constants. 

In the case of ^ = 0.3, the first maximum of the central nucleon density {p^ = 0.30 fm~^) 
is reached at t = 10 fm/c. At later time the system rebounds and oscillates approaching 
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gradually a static configuration with the nucleon density pn — 0.26 fm~^ at the center. 
Since the annihilation is switched off in this calculation, the compressed configuration may 
exist, in principle, infinitely long time. However, due to numerical reasons, the stability is 
destroyed by a gradual test particle escape from a box in the coordinate space, where the 
mean field is computed. Nevertheless, the numerical accuracy is good enough to trace the 
stable system up to at least t = 100 fm/c. 

In the case of ^ = 1, the compression process is much faster than in the case of ^ = 0.3. 
Already at t = 4 fm/c we observe the first maximum of the central nucleon density with 
Pn = 0.48 fm~'^. A smaller value pn ~ 0.34 fm~'^ is reached asymptotically after some 
oscillations. 

In Fig. m we compare time evolution of the nucleon density distribution along the central 
axis z for the light (p^^O) and heavy (p^°^Pb) systems. For p^^O, the bell-like shape of the 
density distribution is reached quite fast. However, in the case of p^°^Pb, we observe a quickly 
growing peak in the center, while peripheral nucleons still do not react on the compression. 
This leads to the delayed shape equilibration via a comphcated compression-decompression 
cycle. 

Fig. [5] shows the time evolution of the central nucleon density for the three systems: 
p^^O, p'^^Ca and p^^^Pb. The case without annihilation is shown by the dotted lines. For 
a comparison, we also present the central density time evolution in the respective ground 
state nuclei without an antiproton inside (dashed lines). We see, that at long enough times 
of the order of several tens fm/c, the static compressed configuration is indeed reached. 
The small oscillations of the central density in the compressed system visible at t > 50 
fm/c are approximately of the same amplitude as the oscillations of the respective ground 



state. Thus, the reason for these small amplitude oscillations is t 



le fermionic ground state 



instability due to the classical treatment of particles [26], l36|, l37|, l38|, |39|. In principle, this 
instability can be removed by either employing the Pauli potential as in Refs. 36|, l37l. l38[or 
by adding a friction force to the Hamiltonian equations of motions for the test particles 39|. 
Such modifications are beyond the scope of the present work. We expect that they would 
not essentially modify the density profiles of the compressed configurations, which are in an 
overall agreement with the previous static Hartree calculations of Refs. 

One can also notice (see lower panels of Fig. [5]) a peculiar feature of the p^^^Pb system: 
the dip in the central nucleon density at t ^ 50 fm/c. This is mostly a consequence of 
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the delayed shape equihbration mentioned above in discussing Fig. HI The especially strong 
density drop for ^ = 0.3 is partly caused by the symmetry loss due to the finite number of 
test particles. The dip is absent in lighter systems, since the shape equilibration for them is 
much faster (10 ^ 20 fm/c, see upper panels in Figs. OandH]). 

Since the compression time is of primary importance, we have also studied the sensitivity 
of our results to the width of the initial antiproton DF (129|) . We have found, that 
for a larger (smaller) width the compression time becomes somewhat longer (shorter). In- 
particular, for the p^^O system at ^ = 1 the time needed to reach the first density maximum 
is 5 fm/c (2.5 fm/c) for (7^ = 2 fm (a^ = 0.5 fm). For the same system at ^ = 0.3 the first 
density maximum is reached at 10 fm/c (3 fm/c) for = 2 fm (a^ = 0.5 fm). The value of 
the maximum density reached in the compression process is practically insensitive to a^. 

We have to admit also, that there is some numerical uncertainty in our calculations due 
to the choice of the width L of the test particle Gaussian and the grid step size. E.g. in the 
calculation with A^,. = = A^ = L = 0.33 fm for the lightest system p^^O the maximum 
and saturation densities are 20 % higher than in calculation with A^^ = Aj^ = A^ = L = 0.5 
fm. Setting the smaller grid step is not feasible for technical reasons. Overall, this numerical 
uncertainty is comparable to the one due to different choices of the scaling factor ^ in our 
calculations. 

B. Post-annihilation dynamics of residual nuclei 

Next, we study the dynamics of a residual nucleus after sudden annihilation of an an- 
tiproton. The annihilation was simulated as described in Sect. HE. For each considered pA 
system and the scaling factor ^, three different annihilation times tann have been chosen: 
They correspond to (i) the early (tann = 0) annihilation from a non-compressed ground- 
state nucleus, (ii) the annihilation at the time moment when the first maximum of the 
central density is reached, and (iii) the late annihilation from an asymptotic compressed 
configuration. 

One can see from Fig. [5l that the annihilation from a non-compressed ground state nucleus 
(thin solid lines) does not lead to significant expansion of the residual nuclear system. The 
central nucleon density stays always below but close to po in this case. In the case of 
annihilation from the compressed configurations (thick solid and dash-dotted lines), we 
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observe that, for the hght systems p^^O and p Ca, the central nucleon density decreases 
sharply after annihilation and reaches values well below po- This is a clear indication of the 
collective expansion of a system from the initially compressed state. On the other hand, 
for the heavy p^^^Ph system, the expansion is not very pronounced at any choice of the 
annihilation time. 

It is interesting, that if the annihilation is switched on at the first density maximum 
(thick solid lines), then after an abrupt falling down the central density stays for some 
time ~ 10 — 20 fm/c close to po before decreasing further. This is explained by an inertial 
compression: After annihilation, the periphery of a residual nucleus still continues to move 
to the center during some time until rebound. In the case of the p^°Ca system, this is 
demonstrated in Figs. M and [TJ where we show the baryon density and the radial collective 
velocity, 

^^rad = r ■ Vcou/r, (34) 

at several times as a function of the radius r. The collective velocity has been determined 
as 

<„u = T"7T°°, « = 1,2,3. (35) 

At r > 4 fm, the evolution of the radial collective velocity field is noticeably influenced by 
the fermionic ground state instability discussed above. Nevertheless, one can still observe 
the inertial compression. Indeed, the radial collective velocity at t = 12 fm/c for C, = 0.3 and 
at t = 6 fm/c for = 1 reveals the fast outward motion at the center, while the peripheral 
nucleons still continue to move to the center. This explains plateaus in the central density 
evolution near t ~ 20 fm/c in Fig. [5l At t 30 fm/c the whole system starts to expand. 
This is reflected in the monotonically increasing radial collective velocity with radius. The 
especially strong rise of frad at large r is due to emission of fast particles. At later times 
t ~ 50^60 fm/c the expansion is replaced by the inward motion of the matter in the central 
zone. However, the fast particles are still continuing to escape from the dense region. We 
expect that in reality the system will break-up into fragments before the inward motion will 
start (see discussion Sect. IV A). 
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IV. OBSERVABLE SIGNALS 



Multifragmentation of residual nuclei 



It is presently well established (see 40|] and refs. therein) that nuclear matter at low 
densities (p < 0.6 po) becomes unstable with respect to small density perturbations, so 
called spinodal instability. However, in order these density perturbations to develop into 
nuclear fragments, the system must stay long enough time ~ 30 fm/c in the spinodal region. 
One can see from Fig. [5l that the light systems p^^O and p^°Ca spend a long time in this 
region. Therefore, the residual nuclear systems can undergo a multifragment break-up, if the 
annihilation happens in the compressed configurations. In other words, the multifragment 
break-up of nuclei after the p-annihilation may serve as a signal of the compression prior the 
annihilation. 

In Table [I we collect the estimated parameters of fragmenting sources for the p^^O and 
p'^^Ca systems. The case of p-annihilation from the state of maximum central density is 
considered here (see thick solid lines in Fig. The sources have been determined by 
selecting nucleons in the space region where the baryon density is larger than p^in = O.Olpo- 
They are characterized by the neutron (A^) and proton (Z) numbers, the collective kinetic 
energy per nucleon {E^°^) and the residual excitation energy per nucleon (E*^^). These 
parameters are defined as 



N= J dhpr,{v) , Z= J rfVpp(r) , (36) 

P>Pmin P>Pmin 

E^^ = j I dh {t'' - ^T^^) , (37) 

P>Pmin 

Els = \ J d'rT''-E,.,.{N,Z)-E^:^ , (38) 



p>Pi 



mm 



where p„ and Pp are the neutron and proton densities, respectively; A = N + Z; and 
Eg.s. {N, Z) is the ground state energy per nucleon of a nucleus with neutron number 
and proton number Z computed within our model. The collective kinetic energy (1371) is 
calculated neglecting pressure effects. Due to the initially isospin-symmetric nuclei and the 
neglect of the isovector mesons in the Lagrangian density ([3]) , we obtained in all cases ~ Z 
in the source. It turned out also, that the Coulomb interaction optionally included in some 
of the calculations does not change this result. Thus, only the charge numbers Z are given 
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in Table HI The time moment for determination of the source parameters has been chosen at 
35 ^40 fm/c, when the central nucleon density is about 1/3 ^ l/2po; i-e. inside the spinodal 
region. The calculated residual excitation energy is typically 6 -h 10 MeV/nucleon that 
corresponds to the temperatures 4 7 MeV. According to the statistical multifragmentation 
model 4l|, in this energy domain the multifragment break-up is the dominating decay 



channel of residual nuclei. The collective kinetic energy is 1.4 -=-2.1 MeV/nucleon in the 
case of p^^O. This is well above the Coulomb energy of the source, which is only about 0.4 
MeV/nucleon for N = Z = 5. Unfortunately, such a source is too small to experience the 
real multifragment break-up, rather a Fermi break-up into small clusters [41']. In the case of 
larger sources, produced in p'^^Ca annihilation the collective kinetic energy is considerably 
smaller, 0.3^0.6 MeV/nucleon, but still significant with respect to the total Coulomb energy 
~ 1.0 MeV/nucleon for N = Z = IQ. Thus, we expect some signs of collective expansion to 
be visible in kinetic energy spectra of produced fragments. 



B. Knock-out nucleon spectra 

Let us now consider other observable effects. Fig. [8] shows the cm. kinetic energy spectra 
of the nucleons emitted from the p^^O, p^°Ca and p^^^Pb systems after the p- annihilation. In 
order to separate emitted nucleons from the bound nucleons of a residual nucleus, we used 
a simple criterion: only those nucleons, both protons and neutrons, were included in the 
spectra which are separated by at least 3 fm from the other test particles of a given parallel 
ensemble at t = 100 fm/c. One can see, that nucleons with the kinetic energy E]^^^ ^ Ep, 
where Ep ~ 35 MeV is the Fermi energy of the nuclear matter at po, are abundantly 
emitted. Such nucleons are knocked-out from the nucleus by the mesons produced after the 

n 

annihilation |14|. 

In Table HTl we list the slope parameters T/v of the nucleon kinetic energy spectra obtained 
by the Maxwell-Boltzmann fit 

dNr, 



dE, 



kin 



(X JEi,i^exp{-Ei,ijTN) (39) 



in the region of i?kin = 200 -^ 500 MeV. We would like to mention, that the authors of Ref. 



14l | report the slope temperature of about 60 MeV for the kinetic energy spectrum of the 



emitted protons in the case of p^°Ca system, which is not so far from our results for the 
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annihilation in the ground-state nucleus at tann = 0. 

We want to emphasize that the kinetic energy spectra of nucleons emitted after the p- 
annihilation from the compressed p^^O, and p'^^Ca. systems are significantly harder than 
the spectra of nucleons from the annihilation at tann = 0. This can be explained by two 
effects. First, the collective expansion of the outer shell will increase the slope temperature, 
typically, by several MeV (see Table [land Figs. [6l [7]). Second, just after the annihilation the 
nucleon potential at the center of a nucleus grows suddenly by ~ 80 4- 300 MeV. This creates 
an additional push for the fast nucleons emitted from the nucleus. Although the hardening 
effect is most pronounced for the lightest system p^^O, it is also quite visible for p^^Ca. For 
the heaviest system p^^^Ph, we observe almost identical high energy tails of the nucleon 
spectra for the different annihilation times. The reason is that the collective expansion is 
practically absent in this system (see Fig. [5]). Also, the yield of fast nucleons is reduced by 
their subsequent rescatterings in the residual nucleus. 

C. Mesonic observables 

The meson production from the p-annihilation inside nucleus is influenced both by the 
mean field via the potentials of the annihilating pair and by the final state interactions 
(FSI), i.e. the two-body collisions and resonance decays. It is instructive to disentangle the 
contributions of the mean field effects from the rescattering and absorption effects. To this 
aim, we have performed additional calculations by subsequently switching off the FSI and 
the mean field. Corresponding results are shown in Figs. [9][T0l and [TTl 

Figure [9] shows pion multiplicity distributions for the p^^O system. In the case of reduced 
antiproton coupling constants (^ = 0.3), the mean field alone does not produce any notice- 
able modification of the pion multiplicity distribution. On the other hand, FSI leads to a 
rather substantial shift toward smaller n^. For the case of ^ = 1, we observe a strong pion 
multiplicity reduction due to smaller ^/s for the pN annihilation, while the FSI effects are 
much weaker. One can also see a quite significant compressional effect for ^ = 1, which is 
only very weak for ^ = 0.3. 

In Fig. [TOl we present charged pion momentum spectra in the cm. frame of the p^^O 
system. FSI strongly modifies these spectra, mostly due to the ttN —>■ A —>■ ttN processes, 
which effectively decelerate pions. As we have already observed earlier in Fig. [HI emitted 
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nucleons gain energy, correspondingly. 

The effect of the baryonic mean field on the pion momentum spectrum is relatively mod- 
erate: we observe some depletion of the high-momentum tail, which is more pronounced in 
the case of ^ = 1. The compressional effect is visible in the reduction of the pion yield in the 
momentum range 0.3 < k < 0.6 GeV/c. This can be understood from Figs. [2] and O in vac- 
uum, the events with 7^,^ = 5 and n,r = 6 contribute substantially to this momentum range, 
while the probability of such events is substantially suppressed in a compressed nucleus due 
to the reduced annihilation phase space 2|. 

The pion momentum spectrum has clearly the two components: the slow pions, which 
have undergone rescatterings via the A-resonance excitation {k less than about 300 MeV/c), 
and the high energy pions, which were emitted from the system almost without secondary 
interactions. The similar result has been obtained in earlier intranuclear cascade calculations 
[l^ . Following Ref. [l^, we have also fitted the low energy part of the pion spectrum 



{E — m^r = 100 -T- 150 MeV, E = ^Jk"^ + ttl^) by the Maxwell-Boltzmann distribution 

'^^^^ eexp{-E/T^) . (40) 



dk 

This fit has produced the following slope temperatures of the charged pion momentum 
spectrum for the p^^O system: ~ 45 MeV and 44 MeV for ^ = 0.3 (T^ ~ 43 MeV and 36 
MeV for ,^ = 1) in the case of early annihilation (tann = 0) and annihilation at the time of 
maximum compression, respectively. The extracted T^r is smaller than the slope temperature 
of 53 MeV of the low energy pions for the p^°Ca system reported in Ref. [ij]. However, in 
the calculation without mean field, we obtain ^ 51 MeV, which is in a good agreement 
with the result of Ref. lUl. Therefore, the difference between our results and those of Ref. 

n 

|14l | is caused by the mean field acting on the annihilating pair in medium. Moreover, we see 
the softening of the pion spectrum in the case of larger antinucleon couplings due to smaller 
y/s of the annihilating pair. 

Figure [TT] presents the distributions of the annihilation events in the total invariant mass 
of produced mesons from the p^^O system. The invariant mass is defined as 

Minv=((PLs)'-PLs)'^' , (41) 



where P^es = S is the sum of four- momenta of the mesons produced in a given annihila- 

i 

tion event. The calculations were done for the case of t^nn = 0. 
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In the absence of FSI, Minv should be equal to the invariant energy ^/s of the annihilating 
pN pair. Indeed, without FSI and without mean field, as expected, we get a quite sharp 
peak at 2mN only slightly smeared out due to the Fermi motion of nucleons and momentum 
spread of the initial antiproton DF ( 129|) . The baryonic mean field leads to the shift of a peak 
position toward smaller Minv and to some broadening of the distribution. The broadening is 
due to the spatial spread of the initial antiproton DF (!29|) . which results in different mean- 
field potentials acting on different annihilating pN test particle pairs. Additionally, the FSI 
leads to a very strong broadening of the invariant mass spectrum due to the deceleration 
and absorption of the annihilation mesons. This is clearly seen in Fig. [TT] for the case when 
the RMF was switched off (dotted line). Nevertheless, the full calculation (solid line) shows 
quite strong softening of the Minv distribution due to the mean-field effects. Obviously, this 
effect is stronger for the case of ^ = 1 as compared with ^ = 0.3. 

Finally, in Fig. [121 we systematically study how the meson invariant mass spectra are 
affected by the nuclear compression effects. The results are shown for the different p-nucleus 
systems. Due to the strong reduction of the nucleon effective mass with the scalar density, 
the meson invariant mass spectra become softer when the annihilation happens in the com- 
pressed configurations, as compared with the annihilation in the normal state at tann = 0. 
The effect is, again, more pronounced for the light systems p^^O and p^°Ca. In the p^^O 
system the shift is almost 500 MeV even in the case of ^ = 0.3. 

We believe that results presented in Figs. ISMTOl and [12] constitute the set of observables 
sensitive to the compressional effects in nuclei induced by an antiproton. 



V. SUMMARY AND OUTLOOK 



We have performed dynamical modeling of possible compression effects in nuclei due 
to the presence of an antiproton. The semiclassical transport GiBUU model 10| incorpo- 
rating the relativistic mean fields for the baryons and antibaryons has been employed in 
calculations. The model reproduces reasonably well the earlier static calculations of bound 
p-nuclear systems 

In this work, we did not consider the stopping process of an incident antiproton in a target 
nucleus. This is a rather complicated problem due to the unknown in-medium cross sections 
of the p-scattering and annihilation. This problem will be addressed in a forthcoming paper. 



22 



Instead, we have assumed, that the antiproton has penetrated to the center of the nucleus, 
stopped there due to an inelastic collision, and then get captured to the lowest energy state. 
Such events should be very rare, with a probability of the order of 10"^ for the central 

u n 

collisions [2| • As proposed in Ref . [21] , formation of bound antiproton-nucleus states can be 
triggered by the emission of fast nucleons, pions or kaons. We have shown, that during the 
time interval of 4 -h 10 fm/c after creation of the initial state the central density of the target 
nucleus grows up to the values of 2 4- 3po depending on somewhat uncertain values of the 
antiproton coupling constants. We expect that the life time of strongly bound antiprotons 
can be long enough to observe this cold compression effect. 

Detailed kinetic simulations of the post-annihilation evolution of residual nuclei have 
been carried out at different assumptions on the annihilation time. It is shown, that the 
p-annihilation in compressed light systems, like p^^O and p'^^Ca, leads to the pronounced 
collective expansion of the residual nucleus, which may result in the multifragment break-up. 
Another clear signature of the nuclear compression is the hardening of the kinetic energy 
spectra of emitted nucleons. On the other hand, the invariant mass distribution of produced 
mesons is shifted to smaller invariant masses due to the in-medium reduction of the nucleon 
effective mass at high scalar density. Similar phenomena are expected also for the case of 
A-nucleus bound states which can be produced via the pp AA reaction on nuclei. 

Another interesting possibility is that the compressed zone of the nucleus might undergo 
a deconfinement phase transition. Then one can expect formation of a quark-antiquark 
droplet with a non-zero baryon number and relatively low temperature 21]. 

Our main assumption in the present study was that the annihilation takes place in the 
central region of a nucleus. The experimental selection of the central annihilation events 
is a difficult problem. No clear trigger condition for such events has been invented so far. 
One suggestion is that the central annihilation events, in average, will be characterised by 
isotropic emission of secondary particles and high fragment multiplicity jisl . [l^ . However, 
further theoretical and experimental efforts are needed to develop a good trigger condition 
for the central annihilation. Despite of these difficulties, we propose to study the above 
predictions in antiproton-nucleus reactions at the future Facility for Antiproton and Ion 
Research (FAIR) at GSI (Darmstadt). 
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TABLE I: Fragmenting source parameters for the different annihilating systems and values of the 
scaling factor ^ of the antiproton coupling constant, tann denotes the annihilation time moment, t 
is the time moment when the source parameters have been determined, pat is the central nucleon 
density. Z, E^^ and , are the charge number, collective kinetic energy per nucleon and residual 
excitation energy per nucleon, respectively (see Eqs. ([36]) . ([37|) and ([38]) ). 
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Including Coulomb interaction in propagation of the test particles. 

TABLE II: Slope temperatures Tn (MeV) for the nucleon kinetic energy spectra (see Fig. [8] and 
Eq. ([39p ). Only the values of T/v for tann = (first number) and for the annihilation at the time 
of the maximum compression (second number) are given. Statistical error is ±2 MeV. 
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FIG. 1: (color online) Pion multiplicity distribution for pp annihilation at rest in vacuum. Data 
points are from Ref. [qI. The dashed line represents the data fit [34] with the Gaussian P{n-,^) = 
exp{ — (n^r— < >)^/2o-,^ }/J2TTa'^^ where < >= 5.01 and cr^ = 1.04. 
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FIG. 2: (color online) Charged pion momentum distribution for pp annihilation at rest in vac- 
uum. The total calculated distribution is shown by the thick solid line. The calculated partial 
contributions from events with various pion numbers are also depicted (see key for notations). The 
calculations are normalized to the number of charged pions per annihilation event. Data from Ref. 
are in arbitrary units and are rescaled to agree with calculations at /c = 0.3 GeV/c. 
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FIG. 3: (color online) Nucleon and antiproton densities (top panels) and potentials (bottom panels) 
vs coordinate z on the axis passing through the center of the p^'^Ca system at selected times 
indicated in the figure. The calculations with the scaling factor ^ = 0.3 (^ = 1) are shown in the 
left (right) panels. Please, notice different scales on vertical axis. 
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FIG. 4: (color online) Nucleon densities along coordinate z at various time moments for p^^O and 
p208pb at ^ = 0.3 and ^ = 1. 
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FIG. 5: (color online) Time dependence of the nucleon density at the center of the p^^O, p'^^Ca 
and p^^^Ph systems for two values of the scaling factor ^ = 0.3 (left panels) and ^ = 1 (right 
panels) of the antiproton coupling constants. The dotted line shows the calculation without p 
annihilation. The thin solid, thick solid and dash-dotted lines show the results with annihilation 
simulated at various times tann indicated in the figure. The dashed lines show the central nucleon 
density evolution for the corresponding ground state nucleus without an antiproton. 
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FIG. 6: (color online) The baryon density (top panel) and the radial collective velocity (bottom 
panel) as functions of the radial distance for the p'^^Ca system computed with ^ = 0.3. The 
annihilation time tann was set to 10 fm/c. 
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FIG. 7: (color online) Same as in Fig. [6l but for ^ = 1 and tann = 4 fm/c. 
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FIG. 8: (color online) Kinetic energy spectra of emitted nucleons in the cm. frame for various pK 
systems and values of the parameter ^. Different histograms correspond to different values of the 
annihilation time tann indicated in the key. 
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FIG. 9: (color online) Pion multiplicity distributions for the p^^O system. The line with full circles 
shows the calculation without mean field and without FSI after the annihilation. The line with 
full boxes shows the result with mean field, but without FSI. Other lines show the full calculation 
at various choices of the annihilation time (shown in the key). The top (bottom) panel presents 
results for ^ = 0.3 (4 = 1). 
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FIG. 10: (color online) Same as in Fig. [9l but for the charged pion momentum distributions in the 
cm. frame of the p^^O system. 
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FIG. 11: (color online) Distribution of the annihilation events on the total invariant mass of emitted 
mesons for the p^^O system. The calculation without mean field and without FSI is shown by the 
line with full circles. The results without mean field but with FSI are shown by the dotted line. 
The line with full boxes shows the result with mean field, but without FSI. The full calculation is 
presented by the solid line. Upper (lower) panel corresponds to ^ = 0.3 (^ = 1). Only calculations 
at tann = are shown. 
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FIG. 12: (color online) Annihilation event distributions on the total invariant mass of emitted 
mesons for various pA systems and values of the parameter ^. Different histograms correspond to 
different values of the annihilation time tann indicated in the key. 



